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Abstract. We have performed a detailed investigation on the world investment networks constructed from 
the Coordinated Portfolio Investment Survey (CPIS) data of the International Monetary Fund, ranging 
from 2001 to 2006. The distributions of degrees and node strengthes are scale- free. The weight distributions 
can be well modeled by the WeibuU distribution. The maximum flow spanning trees of the world investment 
networks possess two universal allometric scaling relations, independent of time and the investment type. 
The topological scaling exponent is 1.17 ± 0.02 and the flow scaling exponent is 1.03 ± 0.01. 
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1 Introduction 

The constituents of a complex system and their interac- 
tions can be characterized by a complex network. The net- 
work perspective has stimulated explosive interests in the 
research of social, infornrational, technological, and bio- 
logical systems, resulting in deeper understanding of com- 
plex systems [T1[21I31I1] ■ As a part of social systems, the 
network properties of many economic and financial sys- 
tems have been studied. This Uterature grows fast and we 
try to present a very brief review below. 

The stocks in a stock market belong to different indus- 
trial sectors. Generally speaking, the prices of stocks in the 
same sector evolve in a correlated manner. If we treat each 
stock as a node and the distance of two stocks based on 
the cross-correlation coefficient as the weight linking the 
two nodes, then the nrarket forms a network. The minimal 
spanning tree extracted from the distance matrix can be 
used to investigate the hierarchical structure of a portfolio 
of stocks, which is usually related to industrial sectors |51 
[51I71IH]- Other topological properties of stock market net- 
works are also studied for different markets [Ql lTOlfTT] . 

Alternatively, rather than considering a portfolio stocks, 
the price time series of a single stock can also be mapped 
into networks, which enables us to investigate the dynam- 
ics of a stock through its network structure. There are 
several methods for this purpose. For a pseudoperiodic 
time series, one can partition it into disjoint cycles ac- 
cording to the local minima or maxima, and each cycle is 
considered a basic node of a network, in which two nodes 
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are deemed connected if the phase space distance or the 
correlation coefficient between the corresponding cycles is 
less than a predetermined threshold [12]. We note that 
a weighted network can also be constructed if the phase 
space distance or the correlation coefficient is treated as 
the weigh of a link. Tlris method for pseudoperiodic time 
series can also be generalized to other time series, which 
has been applied to stock time series [T3]. Other methods 
for network construction from time series are based on 
fluctuation patterns p^llSj. visibility of nodes [16], and 
so on. 

There are also intense interests in the study of world 
trade webs, which describe the international trade rela- 
tions between different economies. Serrano and Bogufia 
have constructed a world trade web utilizing the COM- 
TRADE database of the United Nations Statistics Divi- 
sion and pointed out that the world trade web exhibits 
typical features of complex networks |17j . A fitness model 
[TS] has been proposed to reproduce the topology of world 
trade webs based on a different and more detailed data set 
[TO] , in which the fitness of a node is associated with the 
GDP of the corresponding economy. The fitness model was 
then extended to an evolving and directed description of 
world trade webs [5D]. Furthermore, the interrelation be- 
tween the topology of a world trade web and the GDP 
of countries has been elaborated [21]. The world trade 
web can also be described by weighted networks, in which 
not only the topology but also the trade volume are con- 
sidered [22P23]. A gravity model was used to model the 
weighted world trade web [231113] • The world trade webs 
also show universal allometric scaling [25] . synchroniza- 
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Table 1. The number of nodes A''^, , the number of arcs Na (directed), and the number of edges Ne (undirected) of the world 
investment networks constructed from different kinds of investment data for different years. 



Year 




TP 






ES 






TD 






SD 






LD 




iV„ 


Na 


iVe 


N, 


Na 


N, 


TV,, 


iV, 


N, 


iV„ 


iV, 


iVe 


iV„ 


Na 


iVe 


2001 


195 


3244 


2397 


167 


2261 


1701 


168 


2655 


1959 


120 


971 


778 


168 


2649 


1961 


2002 


194 


3325 


2477 


172 


2315 


1737 


174 


2710 


2042 


136 


1043 


864 


172 


2716 


2044 


2003 


197 


3597 


2649 


177 


2930 


2216 


176 


2910 


2193 


134 


1108 


888 


170 


2893 


2167 


2004 


193 


3791 


2774 


174 


3111 


2344 


173 


3136 


2339 


133 


1156 


930 


165 


2968 


2258 


2005 


195 


4001 


2926 


173 


2764 


2055 


170 


3349 


2482 


130 


1260 


1016 


164 


3117 


2367 


2006 


210 


4483 


3248 


193 


3187 


2353 


182 


3592 


2638 


145 


1353 


1080 


172 


3415 


2528 



tion 
ture 



], community structure [27j . and rich-club struc- 



There are also other economic networks, such as the 
world exchange arrangements web [l^ j the "product space" 
networks [29], the venture capital networks [30ll31| . the 
stock investment networks |32| . and the bank connection 
networks [33|. In this work, we shall study the statistical 
properties of a new economic network named the world 
investment network (WIN). The remainder of the work 
is organized as follow. In Section [51 we briefly describe 
the data sets adopted and construct the world investment 
networks. We investigate the basic statistical properties of 
unweighted networks in Section [3| and weighted networks 
in Section [3| We further study the allometric scaling in 
Section [5] Section [B] concludes. 




2 Constructing world investment networks 



We construct the world investment networks using the 
Coordinated Portfolio Investment Survey (CPIS) data re- 
leased by the International Monetary Fund (IMF). The 
CPIS data are publicly available and can be retrieved from 
the IMF web site. There are five kinds of data including 
total portfolio investment (TP), equity securities (ES), to- 
tal debt securities (TD), long-term debt securities (LD), 
and short-term debt securities (SD). All these data are 
recorded from 2001 to 2006. Therefore, there are 30 tables 
in total. Each table contains the investment information 
among different economies. 

A directed and weighted network can be constructed 
from each table. In each network, the nodes represent 
economies. If economy i invests in economy j, we can draw 
a directed link i —^ j from node i to node j, to which a 
weight identical to the investment volume is assigned. The 
directed and weighted network can be fully expressed by 
a matrix W, where the element Wij stands for the invest- 
ment volume from economy i to economy j. We note that 
Wii = for all economies and Wij = if economy i does not 
invest in j. The matrix W does not need to be symmetric, 
that is, Wij ^ Wji. A schematic diagram is illustrated in 
Figure [1] which is a part of the network constructed from 
the ES data in 2001. When Wij = 0, there is no link from 
i to j in the figure. 



Fig. 1. Diagram illustrating part of the world investment net- 
work obtained from ES in 2001. The numbers 1, 2, 3, 4, 5, and 
6 in the open circles signify Australia, China, Egypt, Luxem- 
bourg, United Kingdom, and United States, respectively. The 
numbers above the edges represent the network weights (in- 
vestment volumes). 



The directed network W can be converted into an 
undirected network G, whose weight is determined by 



gij = gji 



ujt 



(1) 



For convenience, directed and undirected links are called 
arcs and edges, respectively. Let A^^ be the number of 
nodes of network W . Then, the number of arcs of W is 



N,, 



EEiK 



(2) 



where the indicator function I(.t) equals to 1 if a; > and 
otherwise, and the number of edges of G is 



N^ W„ 



N^ Af„ 



N, 



JEEi(5.) = EE% 



(3) 



4 = 1 J = l 



1=1 J=l 



It is obvious that 2Ne ^ Na- Table [T] reports the values of 
A^t,, Na and A'e for all the 30 networks. Roughly speaking, 
for each of the five investment types, the number of nodes 
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Na is almost constant before 2006 and increases sharply 
in 2006, while the number of links {N^ and Na) increases 
gradually from 2001 to 2006 for both directed and undi- 
rected networks, which implies an increasing globalization. 



3 Basic statistical properties of unweighted 
world investment networks 

3.1 Undirected networks 

An undirected and unweighted network A can be extracted 
from an undirected and weighted network G. We note that 
A is the adjacency matrix of G. The element a^ of A can 
be determined as follows, 



= I(ff..) 



(4) 



Speaking differently, a^ = 1 if economy i invests in j or j 
invests in i or both, and Oy ~ otherwise. 

We first investigate the degree distribution of nodes 
for each network. The degree of node i can be calculated 
as follows 

ki = Y^ aij . (5) 

Figure [2] shows the degree distributions for the 30 undi- 
rected networks. For clarity, the data for different types 
of investment are shifted vertically. For the networks con- 
structed from the same investment type, the six degree dis- 
tributions for different years almost collapse onto a single 
curve. For the networks from different investment types, 
the distributions are different from each other. 




Fig. 2. (Color online) Degree distributions of the constructed 
undirected WINs. The data points for LD, SD, TD, and TP 
have been translated vertically by a factor of 10~^, 10~^, 10"'^, 
and 10~* for clarity. The solid lines are the best power-law fits 
to the data. 



An evident feature of the degree distributions is that 
they all exhibit a power-law behavior: 



It means that the world investment networks are scale- 
free. The power-law exponent 7 can be estimated through 
a linear least-squares regression to fit the data in log-log 
coordinates. Table [2] reports the exponents 7 for all the 
networks. It is found that, for each type of networks, the 
exponents 7 slightly decreases along time. The observation 
means that there are more highly connected economies 
recently, also indicating an increase in globalization. 



Table 2. The power-law exponents 7 of the degree distribu- 
tions for undirected and unweighted networks. The numbers in 
the parentheses are the errors magnified by 100. 



Year 


TP 


ES 


TD 


SD 


LD 


2001 


1.03(8) 


1.12(5) 


1.04(9) 


1.46(19) 


1.04(10) 


2002 


1.00(8) 


1.13(4) 


1.04(11) 


1.47(20) 


1.03(12) 


2003 


1.00(9) 


1.01(10) 


1.02(10) 


1.39(16) 


0.99(12) 


2004 


0.95(7) 


0.95(9) 


0.95(10) 


1.32(14) 


0.92(11) 


2005 


0.93(9) 


1.04(8) 


0.91(12) 


1.24(12) 


0.90(13) 


2006 


0.93(8) 


1.06(10) 


0.93(12) 


1.33(14) 


0.91(11) 



The average minimum path length is among the most 
studied quantities of complex networks. Table [3] lists the 
path length of all the networks. We see that the values of 
average minimum path length are quite small. For the TP, 
TD and LD networks, the average minimum path length 
decreases. For the ES networks, the average minimum 
path length does not have a clear trend and reaches a max- 
imum value in 2006. For the SD networks, it is also hard 
to identify an evident trend. It is noteworthy to point out 
that the average minimum path length does not change 
much from one year to another. 



Table 3. The average path lengthes of all the undirected and 
unweighted networks. The numbers in the parentheses are the 
errors magnified by 100. 



Year 


TP 


ES 


TD 


SD 


LD 


2001 


1.89 


1.91 


1.86 


1.99 


1.86 


2002 


1.88 


1.94 


1.90 


2.12 


1.88 


2003 


1.87 


1.86 


1.86 


2.03 


1.83 


2004 


1.84 


1.84 


1.84 


2.01 


1.81 


2005 


1.83 


1.91 


1.80 


1.92 


1.79 


2006 


1.84 


1.95 


1.82 


1.97 


1.79 



p{k) 



(6) 



The clustering coefficient of a node is a measure of the 
cluster structure indicating how much the adjacent ver- 
tices of its adjacent vertices are its adjacent vertices. In 
other words, the clustering coefficient of node is the ra- 
tio of the number of existing edges between its adjacent 
vertices to the number of possible edges between them. 
Table |4] presents the average clustering coefficients for all 
the networks. The average clustering coefficients of net- 
works with the same investment type are almost the same 
for different years. A closer scrutiny shows that, the aver- 
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age clustering coefficient roughly decreases along tinre for 
the ES type and increases for the other four types. 



Table 4. The average clustering coefficients of all the undi- 
rected and unweighted networks. 



Year TP 



ES 



TD SD 



LD 



2001 


0.67 


0.67 


0.66 


0.52 


0.64 


2002 


0.73 


0.66 


0.68 


0.49 


0.66 


2003 


0.71 


0.67 


0.67 


0.50 


0.67 


2004 


0.73 


0.66 


0.66 


0.55 


0.66 


2005 


0.73 


0.67 


0.68 


0.60 


0.67 


2006 


0.74 


0.64 


0.70 


0.60 


0.68 



3.2 Directed networks 

A directed and unweighted network B is the adjacency 
matrix of the corresponding directed and weighted net- 
work W. The element 6y of B can be determined as fol- 
lows, 

bij =I(wy) . (7) 

Speaking differently, bij = 1 if econonry i invests in econ- 
omy j and bij = otherwise. The in-degree k" 
fc°"*, and all-degree 7*^'' are defined by, 



out-degree 






/ , bij , k^ 



ki + ki 



(8) 



Figure[3]illustrates the degree distributions of the directed 
and unweigthcd network B constructed from the ES data 
in 2001. 




Fig. 3. (Color online) Degree distributions of the directed 
and unweighted network B constructed from the ES data in 
2001. The data points for in-degree and out-degree have been 
translated vertically by a factor of 0.1 and 0.01 for clarity. The 
solid lines are the best power-law fits to the data. 



We observe that the probability distributions of the 
degrees are consistent with a power-law behavior: 



where io = in for in-degrees, io = out for out-degrees, 
and io = all for all-degrees. Linear least-squares regres- 
sions give the estimates of the three power-law exponents: 
7^*" = 1.17 ± 0.07 for all-degrees, 7'" = 1.06 ± 0.07 for in- 
degrees, and 0.58 ±0.06 for out-degrees, respectively. The 
exponents 7*^", 7'", and 7°"' for the all the directed and 
unweighted networks are reported in TableO On average, 
the exponents 7^", 7'", and 7°"' decrease along time for 
each investment type. This shows that the investments 
among economies become much denser from year to year, 
also a signal of an increasing globalization from 2001 to 
2006. 



Table 5. The power-law exponents 7^" (top panel), 7'" (mid- 
dle panel) and 7°"' (bottom panel) of the degree distributions 
for all the directed and unweighted networks. The numbers in 
the parentheses are the errors magnified by 100. 



Year 


TP 


ES 


TD 


SD 


LD 


2001 


1.05(6) 


1.17(7) 


1.00(7) 


1.20(14) 


1.00(7) 


2002 


1.02(5) 


1.12(5) 


1.03(7) 


1.19(10) 


0.99(8) 


2003 


1.01(5) 


1.00(6) 


1.00(6) 


1.25(12) 


0.97(8) 


2004 


0.97(3) 


0.98(6) 


0.97(7) 


1.15(6) 


0.94(8) 


2005 


0.97(7) 


1.04(6) 


0.94(8) 


1.12(7) 


0.94(9) 


2006 


0.97(6) 


1.02(7) 


0.96(8) 


1.22(5) 


0.94(7) 


2001 


1.03(8) 


1.06(6) 


0.95(8) 


1.38(7) 


0.95(9) 


2002 


0.83(6) 


1.04(7) 


1.50(31) 


1.42(7) 


0.85(14) 


2003 


0.89(7) 


1.01(9) 


1.05(10) 


1.37(10) 


0.91(11) 


2004 


0.77(7) 


0.91(8) 


0.88(8) 


1.30(7) 


0.83(11) 


2005 


0.77(11) 


1.08(7) 


0.79(13) 


1.20(7) 


0.78(13) 


2006 


0.73(12) 


1.03(8) 


0.75(13) 


1.29(5) 


0.71(12) 


2001 


0.35(15) 


0.58(7) 


0.16(15) 


0.81(17) 


0.17(16) 


2002 


0.42(16) 


0.51(10) 


0.32(11) 


0.75(18) 


0.31(16) 


2003 


0.53(24) 


0.19(14) 


0.19(15) 


0.85(15) 


0.18(14) 


2004 


0.23(17) 


0.18(19) 


0.17(19) 


0.75(11) 


0.19(20) 


2005 


0.31(18) 


0.15(15) 


0.19(16) 


0.91(13) 


0.24(16) 


2006 


0.05(10) 


0.23(8) 


0.16(18) 


0.64(27) 


0.14(27) 



Table 6. The average minimum path lengthes of all the di- 
rected and unweighted networks. 



Year 


TP 


ES 


TD 


SD 


LD 


2001 


4.68 


4.24 


4.16 


3.82 


4.16 


2002 


4.65 


4.48 


4.32 


4.58 


4.30 


2003 


4.50 


4.28 


4.21 


4.26 


4.09 


2004 


4.28 


4.01 


3.94 


4.05 


3.85 


2005 


4.24 


4.09 


3.72 


3.76 


3.72 


2006 


4.51 


4.46 


4.05 


4.13 


3.87 



p(fc-) ~ (fc'°)- 



(9) 



We also report in Table [6] the average minimum path 
lengthes of the constructed directed networks. We find 
that, for each type of investment, the average minimum 
path length gradually decreases from 2001 to 2005, fol- 
lowed by a sharp increase in 2006. This is not surprising 
and can be explained as follows. According to Table [1] 
the number of nodes or economies increases abruptly in 
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2006, which have few links connecting to other nodes. This 
considerably increases the average minimum path length 
of a directed network. The increase of the average min- 
imum path length in 2006 does not mean a weakening 
globalization in 2006. On the contrary, the inclusion of 
more economies in the networks indicates a speedup in 
the globalization. 



4 Basic statistical properties of weighted 
world investment networks 

4.1 Distribution of arc weights 

Figure |4] plots the empirical probability distributions of 
weights for the six directed and weighted networks con- 
structed using the TP data from 2001 to 2006. No clear 
evidence of power laws is observed. We find that the re- 
sults are very similar for the ES, TD, SD and LD data 
(see Figure Al in the Appendix). 




Fig. 4. Empirical probability density functions of weights 
of the directed and weighted networks constructed using the 
TP data from 2001 to 2006. The curves have been shifted for 
clarity. The solid and dot-dashed lines are the best fits with 
the WeibuU and g-exponential distributions. 



We apply the Weibull and the g-exponential distri- 
butions to model the weight distributions [341135]. The 
Weibull probability density Pw{w) can be written as 



Pw{w) = a(3w^ cxp{~aw' ), 



(10) 



and its complementary (cumulative) distribution function 

Cw{w) is 

C^,{w)=exp{-aw'^). (11) 

When (3=1, Pw{w) recovers the exponential distribution. 
When < P < 1, pw{w) is a stretched exponential or sub- 
exponential. When /3 > 1, Pwiw) is a super-exponential. 
The g-exponential probability density Pq{w) is defined by 



and its complementary cumulative distribution function 
Cq{w) is 



Cq{w) = [l + {l-q){-^iw)]-- 



(13) 



Usually, we have q> 1. When (1 — q){—fiw) ^ I, wc ob- 
serve a power-law behavior in the tail Cq{'w) ~ ti;~i/(9~i) 
with a tail exponent of l/(g — 1). 

We adopt the nonlinear least-squares estimator (NLSE) 
to calibrate the weight distributions. The objective func- 
tion in the fitting is ^[Inp(w) — Inp(w)]^ rather than 
^[p(w) — p(w)]^, where p{w) is the empirical data. The 
resultant fits are also illustrated in Figure S) The parame- 
ters of the two models and the corresponding root-mean- 
square values Xw and Xg of the fitting residuals are shown 
in Table [71 It is evident from Figure [H that the Weibull is 
a better model than the q-exponential, which is confirmed 
by the much smaller values of Xq compared with Xw in 
Table [3 We also find that the values of a and (3 remain 
constant from 2001 to 2006. Very similar results are ob- 
served for the ES, TD, SD and LD networks (see Table 
Al and Table A2 in the Appendix). 



Table 7. Estimated parameters {a, 
of fitting residuals {xw and Xg)- 



P, fi and g) and the RMS 



Year 




Weibull 


<?- 


exponential 




a 


/3 


x™ 


M 


q 


Xq 


2001 


0.44 


0.22 


0.0010 


0.05 


2.87 


0.0024 


2002 


0.42 


0.23 


0.0006 


0.06 


2.89 


0.0014 


2003 


0.43 


0.22 


0.0001 


0.06 


3.07 


0.0007 


2004 


0.42 


0.22 


0.0003 


0.04 


2.87 


0.0018 


2005 


0.40 


0.22 


0.0003 


0.04 


2.86 


0.0018 


2006 


0.43 


0.21 


0.0004 


0.05 


3.03 


0.0014 



4.2 Distribution of node strength 

For weighted networks, the node strength is a generaliza- 
tion of the degree, which is defined by 






^ji 



(14) 



The node strength distributions of all the directed and 
weighted networks are illustrated in Figure [5l 

For all the node strength distributions shown in Figure 
[5]) , we see nice power- law behaviors 



p{s) 



(15) 



Pq{w) = /i [1 + (1 - q)i-fiw)] 1-' 



fl2) 



The solid lines are the best fits to the data. The power- 
law exponents 7s are depicted in Table [H These power- 
law exponents are close to each other for each investment 
type. For the TP networks, the exponents exhibit a clear 
decreasing trend from 2001 to 2006. 
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Fig. 5. (Color online) Node strength distributions of the di- 
rected and weighted networks. The curves have been shifted 
vertically for clarity. 

Table 8. The power-law exponents 7s of the node strength 
distributions. The numbers in the parentheses are the errors 
magnified by 100. 



the more efficient is the transportation. Any spanning tree 
can range in principle between two extremes, that is, the 
chain-like trees and the star-like trees. A chain tree has 
one root and one leaf with no branching. For chain-like 
trees, it is easy to show that Ai = i and Ci = i{i + l)/2. 
Asymptotically, we have a power between Ci and Ai with 
the exponent 77 = 2~. For star-like trees of size n, there 
are one root and n — 1 leaves directly connected to the 
root. We have A = C = 1 ior all the leaves and A = n 
and C = 2n — 1 for the root. It follows approximately that 
7] = 1+. Therefore, if the relation (fTSl) holds, the exponent 
should be 1 < 77 < 2. 

For each undirected and unweighted network, a min- 
imal spanning tree can be obtained. The calculated C is 
plotted in Fig. [6] as a function of A for each network. Nice 
power-law relations are observed between C and A. The 
points {A ~ 1,C ^ 1) for the leaves are not shown |42j . 
For each investment type, the data points of the six net- 
works collapse onto a single curve. 



Year 


TP 


ES 


TD 


SD 


LD 


2001 


1.13(4) 


1.12(3) 


1.07(2) 


1.11(4) 


1.08(3) 


2002 


1.13(4) 


1.13(7) 


1.08(2) 


1.13(4) 


1.08(2) 


2003 


1.11(4) 


1.06(2) 


1.06(2) 


1.11(4) 


1.06(2) 


2004 


1.09(3) 


1.11(4) 


1.11(4) 


1.10(3) 


1.04(2) 


2005 


1.08(3) 


1.13(4) 


1.10(4) 


1.07(3) 


1.02(2) 


2006 


1.07(2) 


1.13(4) 


1.10(4) 


1.08(2) 


1.09(4) 



5 Universal allometric scaling laws 

5.1 Unweighted world investment networl<s 

Allometric scaling laws are ubiquitous in networking sys- 
tems such as metabolism of organisms and ecosystems 
river networks, food webs, and so on [36 l [37 ( l38 l l39 p 40 ( l4 H 
|42] . For economic systems, the world trade webs also ex- 
hibits a universal allometric scaling |25| . For a complex 
network, a minimal spanning tree can be extracted. Each 
node of the tree is assigned a number 1. Two values Ai 
and Ci arc defined for each node i in a recursive manner 
as follows: 



A 



E 



A, 



1 



and 



c,= J2 ^^■+^* 



(16) 



(17) 



jeJ(') 



where 3{i) stands for the set of daughter nodes of i and 
Ai = Ci = 1 |40]. The allometric scaling relation is then 
highlighted by the power law relation between Ci and Ai : 



C ^A"^ 



(18) 



We note that not all trees have such allometric scaling, 
such as the classic Cayley trees [43] . 

For spanning trees extracted from transportation net- 
works, the power law exponent jy is a measure of trans- 
portation efficiency [40ll42j . The smaller is the value of 77, 



10 



10 



^10-' 



10" 



10" 




10" 



Fig. 6. (Color online) Allometric scaling relationship between 
S and A. The data points for LD, SD, TD, and TP are trans- 
lated vertically for clarity. The solid lines are the best power- 
law fits to the data. 



Nice power-law behaviors are observed in Figure O A 
linear fit of In C against In A gives the estimate of ?/ for 
each network. The trivial point {A = 1,C = 1) should 
be excluded from the fitting [42] . The resulting exponents 
are listed in Table [9] We find that the exponents arc al- 
most the same for different investment types and different 
years, which means that the power-law allometric scaling 
is universal for the world investment networks. This value 
of rj is comparable to 77 = 1.13 ^ 1.16 for food webs [42], 
but much smaller than r] = 1.3 for world trade webs [55] 
and rj = 1.5 for river networks [40| . 



5.2 Weighted world investment networks 

When studying the world trade webs, Duan proposed a 
framework of flow allometric scaling analysis for weighted 
networks [25], 



fA, 



j6J(i) 



-Wi 



(19) 
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Table 9. The topological scaling exponents 77 for all the undi- 
rected and unweighted networks. The numbers in the paren- 
theses are the errors magnified by 100. 



Year 


TP 


ES 


TD 


SD 


LD 


2001 


1.15(3) 


1.10(4) 


1.16(1) 


1.19(2) 


1.16(2) 


2002 


1.14(3) 


1.15(4) 


1.16(2) 


1.27(3) 


1.15(2) 


2003 


1.18(3) 


1.16(2) 


1.16(2) 


1.20(3) 


1.13(3) 


2004 


1.18(3) 


1.17(3) 


1.17(3) 


1.21(2) 


1.14(3) 


2005 


1.18(3) 


1.14(3) 


1.17(2) 


1.18(3) 


1.13(3) 


2006 


1.19(3) 


1.15(3) 


1.16(2) 


1.19(2) 


1.18(2) 


mean 


1.17(2) 


1.15(2) 


1.16(1) 


1.21(3) 


1.15(2) 



and 



fa = Y. f^i + /^» 



(20) 



ieJ(0 



where J(i) is the set of daughter nodes of i. Wi stands for 
the total weight flowed into node i, fAi is the weighted 
quantity of resources and fCi stands for the weighted 
transferring cost. For the leaves, the carried values of fA 
and fC are identical to the investment volumes. One can 
also expect for certain network that the allometric scaling 
relation exists between fCi and fAf. 



fC ~ M« 



(21) 



in which the exponent ( is called the flow allometric scal- 
ing exponent |25j . 

We adopt this analysis on the maximum-flow spanning 
trees of the undirected and weighted world investment net- 
works. The maximum-flow spanning trees of the invest- 
ment networks can be obtained [25] . Figure [7] plots fC 
with respect to fA in double logarithmic coordinates. For 
each investment type, the data points of the six networks 
collapse onto a single curve, independent of the time. 




Fig. 7. (Color online) Allometric scaling relationship between 
fS and fA. The data points for LD, SD, TD, and TP are 
translated vertically for clarity. The solid lines are the best 
power-law fits to the data. 



Evident power-law behaviors are observed between fC 
and fA for all the weighted networks under investigation. 



The curves seem parallel for different types of networks. 
The flow scaling exponents C, are estimated by the slopes 
of the linear fits of In fC with respect to In fA. The ex- 
ponents C are reported in Table [101 We find that the flow 
allometric scaling exponents are almost the same for dif- 
ferent investment types and different years, which means 
that the power-law allometric scaling is universal for the 
world investment networks. It is interesting to point out 
that this (^ value of the world investment networks is close 
to the flow allometric scaling exponent of the world trade 
webs [25, . 



Table 10. The flow scaling exponents ^ for all the undirected 
and weighted networks. The numbers in the parentheses are 
the errors magnifled by 100. 



Year 


TP 


ES 


TD 


SD 


LD 


2001 


1.02(0) 


1.01(0) 


1.03(0) 


1.07(1) 


1.03(0) 


2002 


1.02(0) 


1.01(0) 


1.02(0) 


1.06(1) 


1.02(0) 


2003 


1.02(0) 


1.03(0) 


1.03(1) 


1.05(1) 


1.02(0) 


2004 


1.02(0) 


1.03(0) 


1.03(0) 


1.04(1) 


1.03(1) 


2005 


1.02(0) 


1.01(0) 


1.03(0) 


1.04(1) 


1.03(1) 


2006 


1.02(0) 


1.01(0) 


1.03(0) 


1.04(1) 


1.02(0) 


mean 


1.02(0) 


1.02(1) 


1.03(0) 


1.05(1) 


1.03(1) 



6 Conclusion 

In this work, we have constructed a new type of economic 
networks based on the Coordinated Portfolio Investment 
Survey data released by the International Monetary Fund, 
which records data from 2001 to 2006. We have studied the 
statistical properties of these world investment networks. 
Our results show that there is an increasing globalization 
in the past few years under investigation. 

The degree distributions are scale-free for all the con- 
structed networks. For the same investment type of data, 
the average path length and average clustering coefficient 
are almost the same for different years with a weak trend. 
When we regard the world investment networks as weighted 
networks, the WeibuU and g-exponential distributions are 
utilized to flt the weight distribution by means of a non- 
linear least-squares estimator. We find that the WeibuU 
model remarkably outperforms the g-exponential model. 
In addition, the node strength distributions arc found to 
exhibit nice power-law behaviors. 

We also investigated the allometric scaling of the mini- 
mal spanning trees and the maximum-flow spanning trees 
of the world investment networks. There are two universal 
allometric scaling exponents characterizing the topological 
structure and the investment pattern of the networks. We 
find that the topological scaling exponent is 77 = 1.17±0.02 
and the flow scaling exponent is C = 1.03 ±0.01. The topo- 
logical scaling exponent is found to be close to that of the 
food webs and smaller than that of the world trade webs, 
while the flow scaling exponent is comparable to that of 
the world trade webs. 
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Fig. Al. Empirical probability density function of weights in the constructed networks for different kinds of investments and 
different years. The solid and dash-dotted lines are the maximum likelihood fits to the WeibuU and g-exponential distributions, 
respectively. 

Table Al. Estimated values of parameters {a, 13) by means of NLSE. The values of Xw have been multiplied by 100. 



Year 




TP 






ES 






TD 






SD 






LD 






a 


P 


Xw 


Q 


P 


Xw 


a 


P 


X™ 


a 


P 


x™ 


a 


P 


Xt^ 


2001 


0.44 


0.22 


0.10 


0.46 


0.23 


0.13 


0.40 


0.23 





0.38 


0.24 


0.01 


0.46 


0.22 


0.17 


2002 


0.42 


0.23 


0.06 


0.46 


0.22 


0.17 


0.42 


0.23 


0.04 


0.42 


0.23 


0.03 


0.42 


0.23 


0.04 


2003 


0.43 


0.22 


0.01 


0.38 


0.24 


0.01 


0.39 


0.23 


0.01 


0.41 


0.23 


0.02 


0.39 


0.23 


0.01 


2004 


0.42 


0.22 


0.03 


0.40 


0.23 





0.40 


0.23 





0.41 


0.23 


0.02 


0.40 


0.23 


0.01 


2005 


0.40 


0.22 


0.03 


0.44 


0.22 


0.10 


0.42 


0.22 





0.44 


0.22 


0.04 


0.43 


0.22 





2006 


0.43 


0.21 


0.04 


0.42 


0.22 


0.12 


0.40 


0.22 


0.01 


0.44 


0.22 


0.04 


0.42 


0.22 


0.01 


Table A2. 


Estimated values of parameters 


(m, <?) by 


means 


of NLSE. The values of 


Xw have been 


muhiplied by 100. 


Year 




TP 






ES 






TD 






SD 






LD 






M 


<? 


Xq 


M 


<? 


Xq 


M 


1 


Xq 


M 


<? 


Xq 


M 


Q 


Xq 


2001 


0.05 


2.87 


0.24 


0.07 


2.80 


0.22 


0.04 


2.74 


0.11 


0.05 


2.91 


0.07 


0.06 


2.84 


0.30 


2002 


0.06 


2.89 


0.14 


0.06 


2.84 


0.30 


0.06 


2.88 


0.14 


0.05 


2.77 


0.17 


0.05 


2.82 


0.15 


2003 


0.06 


3.07 


0.07 


0.05 


2.91 


0.07 


0.05 


2.91 


0.07 


0.05 


2.86 


0.10 


0.05 


2.87 


0.07 


2004 


0.04 


2.87 


0.18 


0.04 


2.74 


0.11 


0.04 


2.73 


0.10 


0.05 


3.03 


0.10 


0.05 


3.08 


0.06 


2005 


0.04 


2.86 


0.18 


0.05 


2.88 


0.24 


0.04 


2.84 


0.10 


0.06 


3.09 


0.13 


0.05 


3.08 


0.07 


2006 


0.05 


3.03 


0.14 


0.05 


2.97 


0.25 


0.04 


2.84 


0.13 


0.06 


3.17 


0.13 


0.05 


3.18 


0.08 



